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ABSTRACT 

A method for measuring the difference between centroids of polarized flux and total 
flux of an astronomical object - polaroastrometry - is proposed. The deviation of the 
centroid of flux corresponding to Stokes parameter Q or U from the centroid of total 
flux multiplied by dimensionless Stokes parameter q or u respectively, was used as a 
signal. The efficiency of the method is demonstrated on the basis of observations made 
in the V band by using an instrument combining features of a two-beam polarimeter 
with a rotating half-wave plate and a speckle interferometer. The polaroastrometric 
signal noise is 60-70 pas rms for a total number of accumulated photoelectrons Ng of 
10® from a 70-cm telescope; this corresponds to a total integration time of 500 sec and 
an object magnitude V = 6 mag. At smaller Ne the noise increases as ~ \.l"j 
while at larger A^e it remains the same owing to imperfection of the half-wave plate. 
For main sequence stars that are unpolarized and polarized by interstellar dust and 
the Mira type variable R Tri the signal was undetectable. For the Mira type variable 
X Cyg the polaroastrometric signal is found to be 310 ± 70 and 300 ± 70 pa& for Stokes 
Q and U respectively; for o Get these values are 490 ± 100 and 1160 ± 100 pas. The 
significant value of the polaroastrometric signal provides evidence of the asymmetry 
of the polarized flux distribution. 

Key words: instrumentation: polarimeters - techniques: polarimetric - stars: AGB 
and post-AGB 


1 INTRODUCTION 

Quite frequently, some fraction of the radiation of an astro¬ 
nomical object becomes polarized as a result of physical con¬ 
ditions and/or its intrinsic asymmetry and can be mea sured 
by means of conventional polarimetry llTinbergenlll996l l. For 
the same reasons the centroid of an object can depend on 
the orientation of the transmitted polarization. 

This effect is akin to the dependence of the centroid on 
the wavelength of observation, the measurement of which 
constitu t es the subject of spectroastrometry, as proposed by 
iBeckersI (Il982l l. The spectroastrometric signal can be mea¬ 
sured on scales smaller than the diffraction limit of the tele¬ 
scope and its analysis allows one to i mpose additi onal con¬ 
straints on th e models of binary stars (lBailevlll998l l. circum- 
stellar discs jP ontoppidan et all I 2 OIII '). and active galactic 
nuclei liGnerucci et al.ll2013l '). 

Similarly, measurement of the deviation of the centroid 
of the polarized flux from the centroid of the total flux - 
polarimetric astrometry, which we call for brevity polaroas¬ 
trometry - may provide essentially new inform ation about 
the object Jjohnson et al.l[20l3 : lsAnovll2013ll . 

It is convenient to characterize partially polarized light 
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by the so called Stokes parameters I, Q, U, and V, which are 
related to the polarization ellipse in the following way: 

1 = 1, 

Q =/pcos2xcos2G 
U = Ip sin 2x cos 2^, 

V = /psin2^, 

where I is the total intensity, p is the fraction of polarized 
light, X is the angle of polarization, and ^ is the angle char¬ 
acterizing elliptical polarization. ^ varies from —7r/4 to 7r/4, 
for linearly polarized light ^ = 0, for circularly polarized 
light 5 = ±7r/4. Eqs. (O show that the parameters Q and 
U define linear polarization and V defines circular one. The 
Stokes parameters are preferable because they have the same 
units of measurement as the intensity. 

The appearance of an astronomical object can be spec¬ 
ified as a dependence of its intensity on the angular coordi¬ 
nate on the sky Oi{a)-, by analogy it is possible to introduce 
the dependencies of other Stokes parameters on a: Oq{ol), 
Ou{ct), and Ov{oi). 

The total Stokes parameters Q, U, and V are integrals 
of the functions Oq(q;), Of/(a), and Ov{oc) over a and are 
measured by means of polarimetry. Usually dimensionless 
















2 B. Safonov 


Stokes parameters are considered: q = Q/I, u = U/I, and 
V = V/I. 

Consider the images of an object obtained by an ideal 
two-beam polarimeter that splits light into horizontally and 
vertically polarized components: 

-Fh(Q:) = 0/(0 ;)-I-Oq(q:), (2) 

Fv(o) = 0/(a)-0Q(a). (3) 

The centroids of these images are: 

^ /o(Oj(o) + OQ(o))da 
** f (Ot(oi) + OQ(a))da 


_ Ja{Oi{cx)-OQ{a))doi 
f{Oi{a)-OQ{oc))doc ■ 

The half-difference of these centroids Aq = (ah — ocv)/2 
can be measured with significant accuracy, as will be proven 
below. The half-difference is connected to the polarization 
properties of the object in a simple way (if one assumes 
1 ): 


^ f aOQ(a)da f aO/(a)da 

Similar values can be defined for the other two Stokes pa¬ 
rameters U and V. The variables ^Q, ^u, and Av that 
constitute the polaroastrometric signal are the basic observ¬ 
ables of polaroastrometry. Let us denote the components of 
these vectors as follows: Aq = (s*,tj), A// = (s^jt^), and 
Av = (sj,tj). 

By knowing the dimensionless Stokes parameters for our 
object, and given that they significantly deviate from zero, it 
is possible to compute the deviations of the centroids of po¬ 
larized flux from the centroid of the total flux, polaricentres 
ocQjCtu, and ay, from the polaroastrometric signal: 


ctQ = Aq/g, au ~ A/j/u, av ~ Av/u. (7) 


The components of polaricentres will be denoted: aq = 
{sq,tq), ctu = {su,tu), and otv = {sv,tv). This derivation is 
correct for Stokes parameters Q,U, and V, but in following 
we will consider only Q and U, which are related to linear 
polarization. 

Let us consider examples of astrophysical objects that 
differ from the observational point of view for the polaroas¬ 
trometry. 


(i) If the distribution of the polarized flux of the object is 
centrally symmetrical (e.g., a main sequence star), or there 
is no polarized light, then the polaricentres coincide with 
the photocentre and the polaroastrometric signal equals zero 
(first column in Fig. [!}. 

(ii) If the object is polarized by interstellar dust, then the 
polaricentres coincide with the photocentr^B (first column in 

Fig.IU). 

(iii) If the source of polarized light is confined within a 
small region (relative to the precision of the method), then 
the polaricentres will coincide with each other, but deviate 
from the photocentre (second column in Fig. [T|. 


^ Strictly speaking, this is correct only in absence of the atmo¬ 
spheric dispersion, see Section 15.21 
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Figure 1. Different types of objects from the point of view 
of polaroastrometry. Shown in rows are Stokes parameters 7, Q, 
U, photocentres (squares) and polaricentres (crosses). Columns 
display the cases of a single star, a single star with a scattering 
object nearby, and a complex object. 


(iv) In other cases the coincidence of polaricentres is pos¬ 
sible only for special and therefore unlikely object geometry 
(third column in Fig.[T}. For complex objects the polaricen¬ 
tres can deviate from the photocentre even when the total 
polarization is close to zero. 

We can conclude that even a qualitative analysis of polaroas¬ 
trometric data - the detection of a deviation of polaricen¬ 
tres from the photocentre or coincidence of polaricentres - 
presents a valuable tool for astronomical object diagnosis. A 
quantitative interpretation of the polaroastrometric signal is 
also possible by using a priori knowledge about an object. 

Here we propose a method of optical polaroastrome¬ 
try and study its performance on the basis of observations 
made with the prototype of the Multimode Fast Camera in¬ 
stalled on the 70-cm telescope in Moscow (see description 
in Section El- The method used to estimate the dimension¬ 
less Stokes parameters q and u from the observational data 
is presented in Section 0 The polaroastrometric signal is 
extracted fr om the same d ata by using differential speckle 
polarimetry (ISafonovll2013l ). Details of the method in an ap¬ 
plication to polaroastrometry and a noise analysis are given 
in Sectional The results in the form of dimensionless Stokes 
parameters and the polaroastrometric signal are given in 
Section [5] for unpolarized stars, stars polarized by interstel¬ 
lar dust and stars presumably possessing a detectable sig¬ 
nal. Section [6] presents conclusions. Some auxiliary results 
are presented in Appendices El and iBl 
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Figure 2. Schematic of the instrument used in the experiment. Numbers are for the following components: (1) pinhole, (2) auxiliary 
lens, forming an image of the pinhole at the first focal plane (with components 1 and 2 being assembled as a single unit which that be 
inserted into the beam manually), (3) linear polarizer, which is inserted into the beam manually, (4) field diaphragm, (5) half-wave plate, 
(6) collimator lens, (7) filter, (8) Wollaston prism, (9) camera lens, and (10) detector. 


2 INSTRUMENT 

For the polaroastrometric experiment we used a prototype 
of the Multimode Fast Camera under design for the 2.5-m 
telescope of the Sternberg Astronomical Institute (SAI). A 
schematic of the instrument is depicted in Fig. [2] The in¬ 
strument combines features of a two-beam polarimeter with 
a rotating half-wave plate and a speckle interferometer. 

We used an Andor iXon 897 electron multiplying CCD 
(EMCCD) camera as the detector. Owing to electron multi¬ 
plication technology, readout noise in this type of detector is 
greatly reduced, at the cost of a twofold increase in photon 
noise (so-called multiplication noise). The frame rate was 
chosen to be 35 frames per second, which is fast enough to 
reduce blurring of the image by atmospheric turbulence. 

The first optical component of the system is a half¬ 
wave plate (HWP), which effectively rotates the polarization 
plane of incoming radiation. This has two goals: measure¬ 
ment of both parameters of linear polarization and calibra¬ 
tion of the instrumental polarization of the optics located 
after the HWP. The HWP (Edmund Optics, EO-49232) con¬ 
sists of several polymer layers selected so that the total phase 
retardance between ordinary and extraordinary beams devi¬ 
ates from TT rad by ^ 0.02 rad inside the V band. A motor¬ 
ized rotation stage Standa 8MRU was used to control the 
position angle of the HWP. 

Two achromatic lenses with a collimated beam between 
them were used as relay optics. These lenses give a magnifi¬ 
cation of « 4.6. A V filter and a Wollaston prism (RIVoptics) 
were installed in the collimated beam. The prism splits the 
beam in two, polarizing it horizontally and vertically. 

A slit diaphragm is installed in the first focal plane; this 
precludes any overlap of the beams split by the Wollaston 
prism and limits the field of view to 13” x 33”. An example 
of an image obtained with the instrument is shown in Fig. El 
The instrument is controlled from a custom command line 
C-I-+ program running under OS GNU/Linux. 

The instrument was installed at the Cassegrain focus 
of the 70-cm telescope AZT-2 located near the SAI build¬ 
ing in Moscow. The equivalent telescope focal length is 
10.5 m; with relay optics of the instrument this is con¬ 
verted to ~ 48.3 m. The angular sampling frequency is 
fd = 3.22 X 10® rad“^, whereas the cutoff frequency of the 
optical system is /c = D/A = 1.27x 10® rad“^. Therefore the 
Nyquist criterion fd > 2/c is fulfilled for this setup, which is 


0 500 1000 1500 2000 2500 

Figure 3. A single frame from Vega, obtained with the proto¬ 
type of the Multimode Fast Camera and an exposure 30 ms. The 
separation of the images is 210 pixels. The brightness scale at the 
bottom below is in counts. 

necessary for correct computation of the Fourier spectrum 
of images. 

We adopted the following observation sequence. The 
detector continuously obtained a series of frames, without 
gaps between them. Meanwhile, the HWP rotated at a con¬ 
stant angular speed. The frame rate and a HWP speed were 
matched so that the HWP position angle 6 changed by 9° 
over one exposure. Because it was not possible to measure 
and record 6 for each frame with this setup, we determined 
it a posteriori. For each series we introduced a linear po¬ 
larizer at a known orientation into the beam for 10-15 per 
cent of the total series duration at its beginning and end. 
During these periods the polarization of incoming radiation 
can be considered as known and it is possible to determine 
corresponding 9. Then 9 can be interpolated in the middle 
of the series, when the linear polarizer was removed from 
the beam. Data obtained using the described technique al¬ 
low the extraction of both polaroastrometry and polarimetry 
information. 

Observations have been conducted in August-October 
2014 in five sessions, with the instrument being removed 
from the telescope between sessions. Throughout a session 
the angular scale w and the position angle </ of the cam¬ 
era were fixed by the equatorial mount of the telescope. We 
measured these parameters for each session for the trans¬ 
formation of measurements from the instrumental refer¬ 
ence system to the equatorial J2000. Because our exper- 
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Table 1 . Angular scale w and the position angle of the camera 
for the five observational sessions, where the position angle is 
being measured eastward from north to the readout register of 
the CCD camera. Dates of the session and object names used for 
the determination of camera parameters are given in the second 
and third column, respectively. 


Session 

number 

Date 

Object 

w, mas/pix 

0,° 

1 

3.08 - 

4.08 

HIP94336 

HIP108917 

68.8 ±0.8 

69.5 ±0.3 

2 

22.08 - 

28.08 

HIP103669 

HIP108917 

68.6 ± 1.7 

34.0 ± 0.6 

3 

11.09 - 

18.09 

HIP94336 

HIP108917 

68.9 ±0.8 

51.6 ±0.4 

4 

19.09 

HIP94336 

HIP108917 

69.3 ±0.8 

322.1 ± 0.3 

5 

02.10 

HIP94336 

HIP108917 

70.4 ±0.7 

51.2 ±0.3 


iment is not very demanding in terms of w and 0 accu¬ 
racy, for their estimation we used binary and double stars 
with well-k nown epheme rides: HIP94336, HIP10 3669, and 
HIP108917 llMasori et al.ll2007l : iFarrington et al.ll2014l l. We 
acquired long-exposure images of these stars and then mea¬ 
sured a separation vector of the stellar images by approxi¬ 
mating one image with the other. By comparing the param¬ 
eters of these vectors with the ephemerides we obtained w 
and (f). The results of our measurements are given in Table 
[T] For each session we also observed at least one polarization 
standard and an unpolarized star. 


3 POLARIMETRY 


reference system I, Qd, Ud, and Vd as follows: 


JLk = Wil + AWkl2)Kk 


X 


I + C[QdCos(4(6'fe -h K(6lfe))) -I- sin(4(6»fc -I- «(6'fc)))] 


-I- C,'Vd5{0k) sin(2(6'fe -I- n{6k))) 


+ VLk, 


( 8 ) 


jRk = Wil - AWk/2)Kk 


X 


I - C[QdCos(4(6'fe -h K{dk))) + [/d sin(4(61fc -|- K,{6k)))] 


- (:'Vd5{ek)sm{2{ek + KiOk))) 


+ VRk- 


(9) 


Here 9k is the mean angle 9 over frame k duration, Kk is the 
mean atmospheric transparency, W is the total transmission 
of the optical system, AWk is the difference in transmission 
for the polarimeter beams, which includes a flat field error. 
The values S{9k) and K{9k) are the parameters of the HWP 
imperfection; their absolute values are relatively small (see 
appendix 1^ for details) v^k and VRk characterize noise that 
is independent for separate frames, e.g., photon noise. 

Coefficients ^ and S are present in Eq. ([8} and (|9} be¬ 
cause 9 as well as its cosine and sine, changes during the 
course of obtaining the frame, and therefore some averag- 
ing occurs. C give s the sense of the polarimeter effectiveness 
llTinbergenlll996l i and equals ~ 1 —(2/3)6y, where 9f is the 
angle through which the HWP passes during the exposure. 
At 9f = 9° the coefficient ^ equals 0.984. 

From Jkk and Juk the Pk values can be computed as 
follows: 


Pk 


Jhk JrU 
Jhk T JRk 


( 10 ) 


The instrument we used is a two-beam polarimeter, whose 
main advantage is that it is insensitive to variations in atmo¬ 
spheric transparency and to stellar scintillation. However, 
it is susceptible to the difference in transmission between 
the two beams of the polarimeter. To mitigate this effect, 
switching of the images corresponding to orthogonal po¬ 
larizations is employed. The switching usually is performed 
with the help of an HWP installed before the polarimeter 
l|Bagnulo et al.llioOfil 'l. In our case, images are switched every 
140 ms. 

One potential drawback of this scheme is its high sen¬ 
sitivity to any imperfection of the HWP. Nevertheless, it is 
possible to estimate this effect by using a generalization of 
the double difference method as described in the following. 
This generalization is closely related to the adopted method¬ 
ology of polaroastrometry (Section |4|. 

Preparatory processing of data is performed as follows. 
First, we average frames with 9 differing by ^ 3°, which is 
small compared to the angle through which the HWP passes 
during one exposure. Then we subtract a bias and a constant 
background from the averaged frames. Fluxes corresponding 
to two images, J^k and jRk, are extracted by using aperture 
photometry, where k is the number of the averaged image. 
The relative positions of the apertures for photometry are 
the same for all frames. 

In Appendix m the fluxes are shown to be connected 
to the Stokes parameters of the object in the instrument 


Substitution of expressions for the fluxes dll) and gives 
P{0) = cos(4i9) + Ud sin(40)) 

oo oo 

+ ujq{n) cos{n9) + oJu{n) sin{n9). (11) 

n=0 n=0 

Here the functions AW, S, and n describing instrumental ef¬ 
fects are decomposed into a Fourier series over 9 and grouped 
into the last two sums. In the following we designate the 
Fourier spectrum of P(0) as Pifg). 

From Eq. (HU one can see the difference between a sig¬ 
nal carrying information about the linear polarization of the 
source and instrumental effects such as variations of the sys¬ 
tem transmission and the HWP imperfection. The signal 
occupies only one frequency in the spectrum Pifg)-, mean¬ 
while, the noise is distributed over all frequencies, including 
the signal frequency. Therefore, the linear polarization can 
be estimated as the amplitude of fourth harmonics of Pifg), 
but it will be inevitably biased by unknown values uju,qi4). 
All that remains is to estimate their amplitudes by interpo¬ 
lation and take these as a measure of the uncertainty of the 
Qd and Ud estimation. 

In Fig. |4]the spectra \Pifg)\^ for three presumably un¬ 
polarized stars are given. As unpolarized stars we have taken 
main sequence stars at dista nces < 66 pc and galact ic lat¬ 
itude > 10°, as was done bv iBerdvugina et al.l (1201 il l. One 
can see that the spectra are quite uniform, and that there 
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Figure 4. Spectrum modulus squared, |F(/e)P normalized by 
the expected level of the photon noise for four objects. Black lines 
are for presumably unpolarized stars; grey one is for weakly po¬ 
larized o Get. The legend gives the total numbers of accumulated 
photoelectrons. 


are no marked frequencies, which demonstrates the possibil¬ 
ity of estimation of the fourth harmonic uncertainty. It is 
notable that \P{fe)'f slowly rises at low frequencies for the 
bright star. This effect can be explained by the fact that the 
photon component of the noise for this star is sufficiently 
low to make the instrumental effects visible. 

The spectrum \P{fe)\^ for o Get (2014 October 2) is 
also presented in Fig. [I] In this case one sees that the fourth 
harmonic signal is significantly higher than the noise, in¬ 
dicating the presence of detectable linear polarization. For 
polarization standard stars the fourth harmonic amplitude 
exceeds the noise level by a factor of 100. 

The advantage of this method of spectral decomposition 
of P (9) over the double difference method llBaenulo et al.l 
|2009^ is that the former enables an adequate estimation of 
the uncertainty of the resulting value caused by instrumental 
effects, especially the HWP imperfection (see Appendix ID). 
For polarimetry such an estimation is important only for the 
brightest stars, but for polaroastrometry the sensitivity to 
the HWP imperfection is signihcantly stronger and a more 
accurate estimation of its effect becomes critical. A similar 
approach to analysing data from polarimeter s with rotating 
HWP s was applied in millimeter waves by ll Johnson et al.l 

l200'it) . 

From Qd and Ud one can estimate the degree of polar¬ 
ization pd and the angle of polarization Xd\ the latter can 
then be transformed to the equatorial coordinate system. 

To verify the method and test its stability we observed 
polarization standards and presumably unpolarized stars in 
every session. The results are given in Table [Din the form of 
p and X- For comparison, values obtained b ^Schmidt et all 
lll992ll are given. For the standards, the degree of polariza¬ 
tion coincides with a precision 3 per cent, and the angle of 
polarization with a precision 1^2°, both of which are ac¬ 
ceptable. The stability of p is 1 per cent, and that of x is 1°; 
the difference between sessions exceeds the errors insignifi¬ 
cantly. 

The small value of the measured polarization for the 


Table 2. Comparison of measured and expected polarization pa¬ 
rameters of standards. The first column lists the number in Henry 
Draper catalogue or Bayer name; the second lists S for observa¬ 
tional session number and the E for value from ISchmidt et al.l 
(Il992ll : the third lists degree of polarization p and the fourth lists 
the polarization angle x- For stars with small polarization, reli¬ 
able estimation of x cannot be made. 


HD / name 

S/E 

p X 10^ 

X.° 

HD7927 

SI 

332.9 ± 0.3 

93.1 ±0.3 

(f) Cas 

S3 

334.5 ± 0.5 

92.7 ±0.4 


S4 

330.3 ± 0.8 

93.1 ±0.3 


E 

329.8 ± 2.5 

91.1 ±0.2 

HD204827 

SI 

547.7 ± 1.8 

59.1 ±0.3 


S2 

548.8 ± 2.8 

60.5 ±0.7 


S3 

547.7 ± 2.7 

59.9 ±0.4 


S4 

549.7 ± 1.7 

60.5 ±0.3 


S5 

554.4 ± 1.7 

61.5 ±0.3 


E 

532.2 ± 1.4 

58.7 ±0.1 

HD184960 

SI 

1.8 ±0.8 


HD221477 

SI 

3.5 ±3.3 


HD210027 

S5 

1.3 ±0.5 



unpolarized stars provides evidence that the instrumental 
polarization of the telescope is small and can be neglected, 
which is expected because the optical design is centrally 
symmetrical. 

The polarimetry results are also given in the form of 
dimensionless Stokes parameters q and u in Table [S] 


4 POLAROASTROMETRY 
4.1 Basics 

To take into account the effect of the HWP the model of 
images in the two-beam polarimeter m and m should be 
extended as follows: 

FDa) = Oi{a) + a{e)OQ{a) + b{e)Ou{a), (12) 


Fnia) = Oi{a) - a{6)OQ{a) - b{e)Ou{a). (13) 

In our case, a{9) = ^ cos(40) and b{9) = ^sin(40), where 6 
is the orientation of the HWP. The polaroastrometric signal 
is 


1 ai + a{d)aQ + b{d)au 

2 1 + a{8)Q + b{e)U 


Oil — a{6)ocQ — b{0)ctu 

I - a{e)Q - b(e)U 


(14) 


By supposing that 1 (which is correct for most 

astronomical objects), it is possible to derive 


A{e) = a{e)AQ + b{e)Au. (is) 

In the following section we show how the values A{9), Aq 
and Au can be estimated from the real data. 

According to a convention accepted in optics, a beam 
with a Stokes vector (7, 7,0, 0) is polarized horizontally, i.e., 
the plane of polarization is oriented along the OX axis. 
In astronomical polarimetry for a such beam the polariza¬ 
tion plane will be oriented towards the north celestial pole 
(NCP). Therefore, for derivations we assume that the OX 
axis is directed towards the NCP. Then, for an object with 
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Sq > 0, tg = 0 the g-polaricentre is displaced towards north; 
for an object with Sq — 0, tq > 0 - towards east, which may 
seem counterintuitive. 


4.2 Extraction of the polaroastrometric signal 

The polaroastrometric signal A can be extracted from the 
observational data in various ways, e.g. by performing a 
simple centroid computation or by approximating one im¬ 
age with another. We solve this task by using a differential 
speckle polarimetry (DSP), a method sim ilar to the double 
difference approach ([Bagnuto e1^1.l [20091 ) . but operating in 
the Fourier domain (ISafonovI 20131 ). 

A principal observable of the DSP method is TZ{f)'- 


nf) = 


^ {FL{f)FUf)) 


(16) 


where / is the spatial frequency, Fl(/) and Ffl(/) are 
Fourier spectra of the images formed by the left and right 
beams of the Wollaston prism, N^r is the averaged number 
of photons forming the right image. The averaging in this 
equation takes place over a finite number of frames. 

The computation of TZ{f) starts with removing bias, 
constant background and equalizing to zero pixels having 
counts < 5 RMS of read-out noise. Then two rectangular 
windows of 100 x 100 pixels each are selected. The left win¬ 
dow is centred on the brightest pixel of the left image; the 
right window is displaced relative to the left by a constant 
vector. The Fourier spectra are evaluated separately for the 
left and right images for many frames and are combined us¬ 
ing Eq. Hi). Only images having 6 differing by < 3° take 
part in the combination. In such a manner a set of TZ{6, f) 
for different 9 values is formed. 

For processing 7i{9, f) data we need a model of its 
formation. As a starting point we use the main result of 
l|Safonovll201.3l ): TZ averaged over the ensemble of samples of 
random atmospheric phase and photon noise is 


7^(/) = (l + A7^(/))7^o(/), (17) 


where TZo is the ratio of visibilities of the object correspond¬ 
ing to left and right beams, and ATZ{f) is derived by the 
differential polarization aberrations in the system. This ex¬ 
pression accounts for three effects; atmospheric optical tur¬ 
bulence, instrumental polarization and photon noise. 

To adequately describe our experiment the following 
instrumental effects additionally need to be taken into ac¬ 
count: (1) deviations of the beams by the Wollaston prism, 

(2) dispersion from the Wollaston prism and the atmosphere, 

(3) differential aberrations of the HWP, and (4) transforma¬ 
tion of the polarization state by the HWP. 

The first two effects are described by multiplying the 
Fourier spectra of images by their respective filters: 

FLif) = (18) 


Fkif) = FR{f)DR{f)e^^^^^^f\ (19) 

where DL{f) and Dii[f) are the filters describing the dis¬ 
persion, and Pfl account for image displacements in the 
left and right beams of the Wollaston prism. Neither effect 
depends on the angle of the HWP as long as atmospheric 


dispersion affects different polarization states equally and 
the Wollaston prism is installed after the HWP. 

After substituting F^if) and Fk{f ) into Eq. (I16II rela¬ 
tion (nzj is modified as follows: 

TZif) = (1 + A77(/))77o(/)D(/)e'""(^-^>, (20) 

where D{f) = Dn(/)/Di{(/) and p = 

By defining the telescope and instrument as two parts 
of an optical system separated by the HWP, it is possible 
to split the impact of differential aberrations into that from 
the telescope, the HWP, and the instrument. 

n{e,f) = TZo{9,f){l + ATZTi9,f) 

+ A77HWp(e, /) -f Ani{f))D{f)F^^<-'’ F^ (21) 

The correctness of t his operation is provided by the additiv¬ 
ity of bias ATZ{f) (ISafonovll2ni3l ). Note that all the effects 
introduced before the HWP and inside it depend on its po¬ 
sition angle 9. 

Let us suppose that some object with TZo 7 ^ 1 is con¬ 
fined to a region much smaller than the diffraction limit 
of the instrument A/D. In this case, decomposing the ex¬ 
ponent in the Fourier transform in a Taylor expansion and 
keeping only t he first term, one obtains (Ijohnson et akllioill : 
ISafono^l2013l) TZo{9) = exp{i7r(4A(0) ■ /)}, moreover, the 
exponent argument is significantly less than unity for all 
frequencies. Substituting A{9) from (I15II . we obtain: 

no{9) = exp{f4<[/,5(0) + fyT{9)]}, (22) 

where 

S{9) = Sq cos(40) + s* sin(46'), (23) 

T{9) = tq cos(40) -I-1* sin(40). (24) 

The values Sq,Su,tq, and t), comprise the polaroastrometric 
signal, as defined above. 

In Appendix m it is shown that the HWP effect 
ATZnwp{9, f) can be represented as 1 -|- A7 ?.hwp(^,/) = 
exp{i(/HWp(6', /)}, while (/hwp < 1. 

As long as the optical system of the telescope is very 
slow (D/64) and possesses central symmetry, for our level 
of pr ecision ATZt can be considered equal to unity dSafonovI 
l2013h . 

By taking all these reasons into account, expression (HU 
can be transformed as follows: 

TZ{9, f) = exp{i4<[/,S(0) + fyT{9)] + i</HWP(e, /)} 

X {1 +Ani{f))D{f)e*^"^'’ F (25) 

Eq. dig) represents the final model of the observational 
data. Note that real observations are not averaged over 
the ensemble and have some noise (which will be discussed 
later): 

n{9k,f) = n{ek,f)+^fkif), m 

where A4(/) is the noise sample. 

Model (I25II allows us to formulate an algorithm for the 
polaroastrometric signal extraction and estimation of its er¬ 
ror from the measurements of TZ{9,f ). The first operation 
is normalizing of TZ{9, f) by its average over the angle 9. 




















Figure 5. Fourier spectra d{0) for o Get and a. Lyr. Black lines 
are for the real part and grey are for the imaginary part. The 
dashed and dash-dotted lines depict level of harmonics RMS es¬ 
timated from the high- and low-frequency parts, respectively. 


Taking into account that the exponent argument in (|25p is 
much less than unity we obtain: 




m /) 

{nd,f))o 


= exp|i47rC[/^S(6») + fyT{9)] + i0HWP,n(6',/)| 


+ V4/), (27) 


where <j)iiWP,n is the normalized value of ^hwp- 

The dependence of phase TZn{0,f) on the frequency 
/ was approximated by least squares with the plane 
4n(^c{8)fx + d{9)fy). To test validity of the model we com¬ 
puted a reduced chi-squared statistic, which is the sum 
of squared deviations of measurements from the model 
weighted by their variances and normalized by the num¬ 
ber of degrees of freedom, which is the difference between 
the number of independent measurements of TZn{9, f) phase 
and number of model parameters. The reduced chi-squared 
turned out to be 1 T 2, which indicates the correspondence 
of the data to the model. 

From Eq. (Hzl) it follows that for the plane parameters: 

c{9) = CS{9) + Gr.i9) + Nx, (28) 


d{9) = CT{9) + Hr,{9) + Ny, (29) 

where functions Gn{9) and Hn{9) characterize the HWP ef¬ 
fect. As a result of the least squares linearity they can be in¬ 
terpreted as parameters of the plane approximating 4>HWP,n- 
Nx and Ny are the components describing the joint effect of 
atmospheric and photon noise. 

After estimation of c{9) and d{9) for different 9 values 
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we decompose them into a Fourier series by this parameter: 
c{9) — s* cos(40) -I- s* sin(40) 

OO OO 

+ ^ gq{n) cos{n9) + ^ ffu(n) cos(n6l), (30) 

n=l n=l 

d{9) = t* cos(40) -I- t* sin(40) 

OO OO 

-I- hq{n) cos{n9) + hu{n) cos{n9). (31) 

n=l n=l 

Here the effect of the HWP imperfection and random noise 
is decomposed in Fourier series Qq, Qu, hq, and From the 
equations one can see that the behaviour of c{9) and d{9) is 
very similar to the dependence of P{9) in the polarimetry 
method described earlier (Section[21 Eq. (IllD i. Therefore we 
analyze them by analogy with Section [3] In both formulae 
the first two components correspond to the signal and the 
sums correspond to the noise input. The signal is estimated 
as a complex amplitude of fourth harmonics and its uncer¬ 
tainty is estimated by interpolation of the noise spectrum 
between adjacent frequencies. 

In Fig. [5] the Fourier spectra d{9) for o Get and a Lyr 
are presented; total number of accumulated photoelectrons 
is 1.5 X 10® and 3.8 x 10^^, respectively. In the spectrum for 
a Lyr, it is easy to distinguish two regions - of low- and 
high-frequency, divided by approximately the fifth harmon¬ 
ics. The high-frequency part results from photon noise and 
other uncorrelated noises, whereas the low-frequency part 
can be explained by the HWP imperfections (see Section 
S3] and Appendix IH). For the low-frequency noise domi¬ 
nance, to estimate the fourth harmonic amplitude uncer¬ 
tainty one should interpolate the noise from adjacent fre¬ 
quencies (dash-dotted line in Fig.SJ. For faint objects, such 
as o Get, it is possible to use a mean high-frequency noise 
amplitude for this purpose. Taking this into account makes 
clear that a Lyr does not show a distinct polaroastrometric 
signal whereas the o Get does. In the next subsection we an¬ 
alyze briefly the properties of the noise components of the 
spectra c{9) and d{9). 

The obtained values and and their errors 

are converted to the equatorial coordinate system using re¬ 
lations HH) from Appendix m The results and their discus¬ 
sion are given in Section [5] 


4.3 Polaroastrometric signal noise 

As a first step we identify the dependence of noise estima¬ 
tions from the high- and low- frequency parts of the spectra 
c{9) and d{9) on the object brightness. The low- and high- 
frequency errors are denoted as cthf and ctlf, respectively. 
For analysis we use the data for the relatively bright star 
5 Gyg (see Section 15.11) . Additional Poisson noise is intro¬ 
duced into the raw data, effectively dimming the object by 
2, 4, and 6 magnitudes. After this, the processing described 
in the previous section is applied and cthf and ctlf are es¬ 
timated. The results are given in Fig. [S] for X and Y coor¬ 
dinates separately. In this subsection it is assumed that the 
X coordinate is measured along the GGD rows. 

As one can see, cthf depends on the photoelectrons 
number Ne as pix (dash-dotted line in Fig. ©, 

which illustrates the fact that the high-frequency part of 
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A unpolarized HF X 

G-O 5Cyg HFX 

A unpolarized HF Y 

G- G 5 Cyg HF Y 

A unpolarized LF X 

5 Cyg LF X 

A unpolarized LF Y 

5 Cyg LF Y 

V polarized HF X 

X 5 Cyg splitted HF X 

V polarized HF Y 

X 5 Cyg splitted HF Y 

▼ polarized LF X 

* 5 Cyg splitted LF X 

▼ polarized LF Y 

* 5 Cyg splitted LF Y 

O MiraHFX 


O MiraHFY 


♦ MiraLFX 


♦ Mira LF 




Figure 6. Dependence of the polaroastrometric signal A uncer¬ 
tainty on the total number of accumulated photoelectrons. The 
black symbols are for the X component, the grey ones are for the 
Y component, empty symbols are for the high-frequency error, 
and filled ones for the low-frequency error (see text). The dash- 
dotted line is for the photon noise level of astrometry 25pix/v/]Ve. 


the spectra c{9) and d{6) is produced by photon noise. It is 
instructive to compare this wi th the photon lim it of astrom¬ 
etry noise, which is jLindeerenlIl978l ). where /3 is 

the full width at half maximum (FWHM) of the image. As 
long as light is split into two halves in our method, and then 
the half-difference between image centroids is taken, for the 
considered signal the photon limit of astrometry noise is also 
/?/%/Ae. Hence the high-frequency noise of the polaroastro¬ 
metric signal determination corresponds to the photon limit 
of astrometry noise at /? = 25 pix or j3 = 1.7", which is 
comparable to the FWHM of a long-exposure image. 

Photon noise dominates at < 10®; for our telescope 
and instrument this corresponds to a total integration of 
500 sec from an object 1/ ss 6 mag. For larger Ne. the low- 
frequency part ctlf gives a larger input. The latter varies 
from 0.6 to 1.5 mpix and does not show a clear dependence 
on Ae. 

Let us consider the correlation properties of the errors, 
for which the series for 5 Cyg is split into five equal parts, 
with (thf and ctlf being estimated for each separately (see 
Fig. ID). From the figure it is evident that the high-frequency 
noise is uncorrelated on timescales comparable to the se¬ 
ries duration and can be reduced by increasing the accumu¬ 
lation time. Meanwhile, for the low-frequency noise this is 


not correct. Therefore we can conclude that a sample of the 
low-frequency error within a single series is approximately 
constant. 

As described in Appendix El the low-frequency error 
can be explained by the imperfection of the HWP given 
that the centre of the beam passing through it deviates from 
its centre by 1 mm. Such a decentring value is to be ex¬ 
pected because the mechanical structure of the instrument 
is susceptible to flexures. This also explains why the low- 
frequency error is highly correlated for a single series of some 
object but changes after pointing to another object. 

In Fig. 1 ^ cthf and ctlf for some objects considered in 
Section [5] are given. As one can see polarized stars and Mira 
type variables exhibit error properties similar to those of 
unpolarized stars, except for sporadic cases. 


5 RESULTS 

To test the method efficiency we observed objects of several 
types: presumably unpolarized stars (type N), stars polar¬ 
ized by interstellar dust (type P), and Mira type pulsating 
variables (type M). The observation method and data pro¬ 
cessing were identical for all objects. Observation circum¬ 
stances and processing results (polaroastrometric signal and 
Stokes parameters) are given in Tabled 

5.1 Unpolarized stars 

Main sequence stars presumably have a very small polaroas¬ 
trometric signal owing to their symmetry and low polariza¬ 
tion even close to their limb. We observed three single stars 
(HD184960, HD221477, and HD9562), satisfying the follow¬ 
ing criteria: spectral class from A to M, luminosity class V, 
distance to the Sun < 66 pc and galactic latitude > 10°. 
HD9562 exhibits a small total polarization which could be 
caused by interstellar dust (an effect som etimes detectable 
even for nearby stars (iBailev et al.ll20i^ ). 

Three additional observed stars satisfy only some of the 
listed criteria and have peculiarities but nevertheless they 
can be classified as unpolarized stars. HD180161 is a BY 
Dra type variable of spectral class G8V and has an angular 
diameter of 330 pas. Brightness variations for BY Dra vari¬ 
ables are caused by starspots occupying a large part of the 
stellar surface and by rotation of the star. Polarization at 
the limb of its disc ca n be crudely assumed to be the same 
as solar, i.e., ~ 10“® (iFluri &: Stenfldll999h . Therefore for 
the most favourable geometry the expected total polaroas¬ 
trometric signal is 0.3 pas, which is much smaller than the 
error of our measurements. 

HD210027 = L Peg is a spectroscopic binary system. The 
angular size of the semi-major axis of its orbit is 10.3 mas, 
its period is 10.2'*, a nd the spectral class es of its components 
are F5V and G8V llKonacki et al.ll2010l) . One can expect a 
small total polarization resulting from its ellipticity but we 
did not detect it. Its polaroastrometric signal is likely to be 
negligible given the object’s symmetry. 

The main component of the spectroscopic binary 
HD186882 = 5 Cyg is brighter than its secondary by 3.4 mag; 
the former has a spectral class of B9.5IV and falls into the 
class of rapidly rotating stars. The expected angular diame¬ 
ter is 550 pas, and the expected degree of polarization at the 
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limb is < 5 per cent llBochkarev et al1ll985h . Therefore the 
total expected polaroastrometric signal for this star is < 25 
pa s, which is undetec table with our instrument. According 
to lBailev et al.l (l20ld i the total polarization for HD 186882 
is (1.08 ±0.02) X 10"'‘. 

For the six mentioned stars the expected polaroastro¬ 
metric signal is much less than the errors of our measure¬ 
ments. This was proven in the experiment, see Table|3] This 
allows us to argue that our measurements do not have any 
bias at least down to 100 pas. A negligible value of bias of 
this kind is expected if one takes into account the central 
symmetry of the optical scheme. 

5.2 Stars polarized by interstellar dust 

An undetectable polaroastrometric signal is expected for 
stars analogous to type N but polarized by interstellar dust 
(type P), because the latter changes the polarization of each 
point of a star image identically. However, a small signal 
can be generated for observations made through the Earth’s 
atmosphere owing to the joint effect of dependence of the 
degree of polarization on the wavelength and atmospheric 
dispersion. 

Indeed, for some polarization standards variation of the 
degree of polarization across t he V band reaches 0 .2 ^ 0.3 
per cent in absolute measure (ISchmidt et al.iri992l l. Hence 
the polarized flux has an effective wavelength differing from 
the effective wavelength of the total flux by AA. For com¬ 
puting the_effective_wavelength we used the stellar spec¬ 
tra of iGunn fc Strvkeri lll983ll . the sensitivity curve of the 
detector, the transmission curve of the filter, and the de¬ 
pendence of degree of the p olarization on the wa velength 
for the considered stars from ISchmidt et all (Il992ll . Differ¬ 
ences in the effective wavelengths are AA = —1.027 nm and 
AA = —0.480 nm for HD204827 and HD7927, respectively. 
These stars also display a dependence of the a ngle of polar¬ 
ization on the wavelength (ISchmidt et al.lll992l l but this can 
be neglected since its maximum amplitude over the V band 
is 0.8°. 

Atmospheric dispersion displaces the polaricentres and 
this displacement is the same for polaricentres q and u, 
as long as AA is the same. The total displacement is 
= \/s't ± tu = £ ■ ■ tana, where e is the an¬ 

gular atmospheric dispersion (e = 2.12 x 10“® rad/nm for 
place and time of observations) and a is zenith distance. 
The amplitude of polaroastrometric signal can be estimated 
according to its definition. At a = 45° for HD204827 and 
HD7927 the amplitudes are 254 and 69 pas, respectively. In 
our experiment such a small signal cannot be reliably mea¬ 
sured. 

5.3 Mira type variables 

Mira type variables are pulsating stars at late stages of evo¬ 
lution. These objects are expected to demonstrate detectable 
deviations of their polaricentres from the photocentre for 
two reasons. Firstly, a significant fraction of the visible light 
radiation is being scattered by the du st envelope in the close 
vicinity of a star dNorris et al.ll20i^ . which leads to polar¬ 
ization of radiation. Secondly, interferometric observations 
sho w that some Mira variables display significa nt asymme¬ 
try llThompson et al.|l2002l : [Ragland et ahlbood l. 


We considered three Mira variables, suitable for obser¬ 
vations in the northern hemisphere from August to October 
of 2014: y Cyg, R Tri, and o Get. The measured total po¬ 
larization and polaroastrometric signal are given in Table [S] 
For R Tri we did not detect a significant signal; however, 
the total polarization for this star is largest in comparison 
with the other observed Mirae. This can be explained by 
the axial symme t ry of this object detected previously by 
iThompson et all (l2002lf using the Palomar Testbed Inter¬ 
ferom eter in near infrared and confirmed by I Ragland et all 
(l2006tl using the Infrared Optical Telescope Array in the 
same region of the spectrum. 

X Cyg and o Get show significant deviation of val¬ 
ues Sq,sl^,tq, and from zero, indicating the asymme- 
try of image, t he sa me conclusion had b een reached by 
[RaglancLgLalJ 1200^ f or y Cyg and by iKarovska et all 
(1 1997^ : [chandler et all (1200^ 1 for o Get (in the visible and 
mid-infrared). Moreover, for these objects the polaricentres q 
and u clearly do not coincide: sj/q ^ sl^/u and tj/q yf r!^/u. 
Therefore the polaroastrometric signal cannot be explained 
by a single polarized source. 

For o Get the situation is complicated by the presence 
of the nearby component M ira B, an accreting white dwarf 
(ISokoloski Sz BildsterJ 120101 1 . which in principle can be po¬ 
larized. On 2014 October 3, 21:55 UT the components of the 
u-polaricentre were measured to be Su = —570±80 mas and 
tu = 500 ± 80 mas; in the same reference frame the position 
of Mi ra B was Sc = —70 mas and G = 480 mas llPrieur et all 
I 2 OO 2 II . Therefore the u-polaricentre was shifted in roughly 
the same direction as Mira B, which can be interpreted as 
the effect of the polarized emission of Mira B seen against 
the background of the “random” polaroastrometric signal 
generated by the Mira A atmosphere. A more definitive con¬ 
clusion on polarimetric properties of Mira B in the visible 
requires more observations conducted in different phases of 
the pulsations. 

To test the methodology we repeated observations of 
y Cyg and o Get in nearby epochs; corresponding measure¬ 
ments agree well (see Table[3]). In case of o Get, between the 
second and third measurements the instrument was rotated 
counterclockwise with respect to the telescope by 57.5°. As 
one can see Sq,Su,tJ, and transformed to the equatorial 
coordinate system, remained the same after rotation. For 
a larger difference in epochs - 55.O'* (between second and 
third measurements of y Cyg), the polaroastrometric signal 
changes significantly. It is worth noting that the brightness 
of y Cyg for the same period decreased from V = 7.5 mag 
to V = 8.5 mag (AAVSO data). 


6 CONCLUSIONS 

The goals of the present work are to develop a methodology 
and to study the precision limits of polaroastrometry. Anal¬ 
ysis has been performed on the basis of observational data 
obtained with the help of the Multimode Fast Camera pro¬ 
totype. By analysing the dependence of the half-difference 
of the photocentres of the images split by a Wollaston prism 
on the position angle 9 of the HWP we obtained estima¬ 
tions of the polaroastrometric signal, i.e., displacements of 
the centroids of the polarized flux Q or U from the pho- 
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Table 3. Circumstances of the observations, the polaroastrometric signal and the dimensionless Stokes parameters for some observed 
objects. Column 1 lists the number in Henry Draper catalogue or Bayer name or Argelander name; Column 2 lists the universal time 
MM-DD HH:MM for 2014; column 3 lists the session number; column 4 lists the mean zenith distance during the series; column 5 lists the 
mean long-exposure FWHM during the series; column 6 lists the total number of accumulated photoelectrons divided by 10®; columns 
7-10 lists the polaroastrometric signal in J2000 coordinates; and columns 11 and 12 list the dimensionless Stokes parameters in J2000 
coordinates. 


HD/name 

1 

UT 

2 

S 

3 

4 

/3," 

5 

N' 

6 

S*, fl&S 

7 

t*, /i.as 

8 

s*, /ras 

9 

10 

q X 10"‘ 

11 

u X 10'‘ 

12 

186882 

08-03 

19:53 

1 

11.7 

2.2 

69 

unpolarized stars 
+20ib30 +10ib50 

±40±30 

±0±50 

0.8T0.5 

0.8T0.5 

184960 

08-03 

20:12 

1 

4.4 

2.2 

5.6 

±130±90 

-70±100 

-0±90 

-90±100 

-1.0T0.8 

-0.5T0.8 

180161 

08-04 

18:42 

1 

9.3 

1.8 

3.1 

±50±70 

±70±80 

±120±70 

-10±80 

-4.7T0.9 

3.4T0.9 

221477 

08-04 

22:45 

1 

24.8 

2.7 

1.8 

-40±180 

±0±180 

±100±180 

±90±180 

0.3T1.8 

-3.9T1.8 

210027 

10-02 

19:51 

5 

32.2 

2.0 

97 

±110±80 

±10±90 

-50±80 

-70±90 

0.3T0.4 

-1.3T0.4 

210027 

10-02 

20:11 

5 

33.6 

2.1 

96 

±90±80 

±20±90 

±10±80 

-100±90 

-0.6T0.3 

0.3T0.3 

9562 

10-02 

21:16 

5 

63.9 

3.9 

18 

±70±100 

±10±100 

±210±100 

-40±100 

3.5T0.5 

-1.4T0.5 

204827 

08-04 

20:35 

1 

12.5 

2.0 

2.3 

stars polarized by interstellar dust 
+90ibl20 +110ibl20 +20ibl20 

±50±120 

-258.2T5.3 

483.0T3.2 

204827 

09-19 

19:00 

4 

3.2 

2.1 

1.8 

±70±180 

-20±160 

-140T180 

±210±160 

-282.7T5.6 

471.3T3.7 

7927 

09-18 

16:45 

3 

46.9 

2.0 

73 

±10±70 

-60±80 

-30±70 

-150±80 

-333.1T0.7 

-31.1T4.3 

7927 

09-19 

17:44 

4 

40.1 

2.2 

71 

-100±70 

±200±60 

-100±70 

±120±60 

-328.4T1.2 

-35.6T4.0 

X Cyg 

08-03 

21:24 

1 

24.5 

2.6 

9.5 

Mira type variables 
+230ib90 +210ib80 

-50±90 

-240±80 

-22.5T0.8 

-28.9T0.7 

X Cyg 

08-04 

19:16 

1 

25.9 

2.3 

9.6 

±220±70 

±220±50 

-0±70 

-320±50 

-23.9T0.6 

-31.5T0.6 

X Cyg 

09-18 

18:05 

3 

23.6 

2.0 

6.5 

±500±100 

±430±110 

-60±100 

±50±110 

-23.1T1.1 

-40.6T1.0 

R Tri 

09-18 

21:45 

3 

30.6 

1.8 

31 

-50±40 

-60±40 

-30±40 

-180±40 

20.5T1.0 

-74.1T0.5 

R Tri 

09-20 

21:37 

4 

31.1 

2.4 

66 

±20±70 

±80±70 

±30±70 

±30±70 

21.0T1.1 

-75.7T0.7 

0 Get 

09-20 

23:01 

4 

60.0 

3.2 

15 

-150T110 

-470T100 

±890±110 

-790T100 

-1.5T0.7 

-9.5T0.7 

0 Get 

10-03 

21:55 

5 

60.2 

2.9 

8.8 

±10±110 

-430T120 

±850±110 

-750T120 

0.9T0.8 

-15.0T0.8 

0 Get 

10-03 

22:37 

5' 

58.8 

3.1 

7.2 

-130T350 

-440±240 ■ 

+1050T350 

-740T240 

1.4T1.4 

-17.0T1.4 


tocentre of the total flux, multiplied by the corresponding 
dimensionless Stokes parameter, q or u. 

We performed an analysis of the polaroastrometric sig¬ 
nal error dependence on circumstances of an experiment. It 
was found that for a total number of accumulated photo¬ 
electrons Ne < 10® the photon noise dominates and can be 
approximated crudely as X.l"jsfW^. At larger Ne. the input 
of imperfections of the given HWP starts to prevail, having 
an approximately constant level of 50-100 fias. This effect 
can be reduced by an order of magnitude if the HWP rota¬ 
tion and beam centre are aligned with a precision of 0.1 mm 
(see Appendix [All. 

Observations of unpolarized main sequence stars, which 
presumably do not have any detectable polaroastrometric 
signal, showed that the measurements are unbiased at least 
at the level « 100 fj,as. Main sequence stars polarized by in¬ 
terstellar dust also do not display a detectable signal demon¬ 
strating the absence of cross-talk between a strong polariza¬ 
tion and the polaroastrometric signal. 

For method verification purposes we selected three Mira 
type variables x Cyg, o Get, and R Tri because they are ex¬ 
pected to have detectable signal. For x Cyg the measured po¬ 
laroastrometric signal is sj = -1-220±70 jras, = ±220±50 
pas, Su = 0 ± 70 pas, and t* = —320 ± 50 pas in the 
J2000 coordinates for 2014 August 4, 19:16 UT. For o Get 
these values are sj = ±10 ± 110 pas, t* = —430 ± 120 pas, 
Su = ±850±110 pas, and t* = —750±120 pas for 2014 Octo¬ 
ber 3, 21:55 UT. A significant polaroastrometric signal gives 
evidence of asymmetry of the polarized visible flux. Asym¬ 
metry wa s detected for x Gyg and o Get earlier in the in - 
frared by I Ragland et al.l J 2 OO 6 II and iGhandler et al.l ll2007h , 


and f or the latter in the visible also by iKarovska et al.l 
lll997ll . For R Tri the signal was not found. 

As a test of the repeatability of measurements we ob¬ 
served these objects at different, but close epochs and at 
different orientations of the instrument relative to the tele¬ 
scope. We did not detect a significant difference in measure¬ 
ments. Also as a by-product we obtained the polarimetry of 
the mentioned objects with a precision of (1 ± 5) x 10“"^. 

Mira type variables are suitable as test objects because 
apparently they have the greatest polaroastrometric sig¬ 
nal among all celestial objects. Though the obtained data 
demonstrate the effectiveness of the method, their astro- 
physical interpretation is obstructed by the complexity of 
these objects. 

Possible applications of polaroastrometry using a large 
telescope to polarized objects described by a simple geomet¬ 
rical model, e.g., some types of circumstellar environments, 
active galactic nuclei, including Sgr A*, and gravitational 
microlensing are very promising. 
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APPENDIX A: INFLUENCE OF HALF-WAVE 
PLATE IMPERFECTION 

Let us suppose that the HWP imperfection can be mod¬ 
eled as fluctuations of a phase retardance <5 and fast axis 
orientation 6 over its area. To derive these fluctuations for 
the given HWP we illuminated it through a linear polarizer, 
then passed the resulting light through another linear po¬ 
larizer, which we could rotate. With the help of a relay lens 
the HWP surface was imaged on the detector. For each point 


on the HWP, values 5 and 6 were extracted from the depen¬ 
dence of the flux on the second linear polarizer orientation. 
The results are given in Fig. \M\ 

As one can see from figure, fluctuations of 6 reach 2°, 
and become more pronounced closer to the border of the 
plate. For 5, the amplitude of the fluctuation is of the order 
of 3 T 5° in a central part of the HWP. On the basis of these 
measurements and with the help of the model developed by 
I Safonov! (120131 ') we estimated 7?.hwp(S, /) (for the definition 
see Section S21) for different beam parameters. Differential 
aberrations produced by the HWP are purely phase aberra¬ 
tions, and the same is true for 7?.hwp(6', /). The RMS of the 
77 hwp(S, /) phase in area |/| < fc is 0.02 rad, in |/| < 0.2fc 
it is 0.002 rad. 

By approximating the phase of 77.hwp(0, /) with a plane 
we evaluated functions G{6) and H{9) (see Section [4.211 and 
then computed their Fourier spectra. These are given in Fig. 
IA2l for a beam diameter 2 mm and different decentring pa¬ 
rameter in comparison with the analogous spectrum mea¬ 
sured in the polaroastrometric experiment for a Lyr. It is 
evident that the observed behaviour of the spectra can be 
explained by the described HWP imperfection at decentring 
~ 1 mm. 

The useful signal corresponds to a complex amplitude 
of fourth harmonics of the Fourier spectra c(9) and d{0). 
The amplitude of the noise of this harmonic is given in Fig. 
IA3I as a function of beam diameter and decentring. As one 
can see, reducing the noise to the level of 0.1 mpix requires 
aligning the beam and the centre of the HWP rotation with 
an accuracy of 0.1 mm, it is also worthwhile to use beams 
of smaller diameter (< 2 mm). 

To derive the HWP imperfection influence on the po¬ 
larimetry let us take its Mueller matrix in the following form 
(where we have kept only small first-order values and inte¬ 
grated components over the beam): 


Mhwp = 


/ 

1 

0 


0 

cos[4(0 -1- k(0))] 


0 

sin [4(0 -1- «:(0))] 

V 

0 

—(5(0) sin[2(0 + k{9 


0 


0 


sin[4(0 + k(S))] S{9) sin[2(0 J- k{9))] 

— cos [4(0 -I- «;(0))] —S{9) cos[2(0 -)- «;(0))] 

J(0)cos[2(0-f k(0))] 1 


(Al) 


In accordance with the specification, the Wollaston 
prism can be considered ideal in our experiment. The first 
rows of the Mueller matrix for left and right beams of the 
prism are 


Mwl = (l 1 0 O) , 

(A2) 

Mwr = (l 0 O) 

(A3) 

(we are interested in first rows only because we can measure 
the intensity only). 

The total Stokes vectors of images for the left and right 
beams can be expressed through the matrices 

S'out.L = MwLMuWP'S'in, 

(A4) 

■Sout.R = MwRMnWpS’in, 

(A5) 
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Figure Al. Fluctuations of the basic parameters of the HWP over its area. The left panel is for the fast axis orientation; the right 
panel is for the contrast at crossed linear polarizers where S' is the deviation of the retardance from tt. 



Figure A2. Spectrum of black and grey lines show the 

real and imaginary parts, respectively. Bold lines are for measure¬ 
ments of a Lyr, and thin lines are for what is expected when the 
HWP is used. Decentring values 0, 1, and 2 mm are shown by 
solid, dashed, and dash-dotted lines, respectively. 



Figure A3. Dependence of the fourth harmonic amplitude of 
spectrum Hn{0) on the decentring parameter for three values of 
beam diameter on the HWP: 2 (black solid line), 3 (dark grey 
dashed line), and 4 mm (light grey dash-dotted line). 


where Sin = (/in, Qin, Uin, Vin) is the Stokes vector of incom¬ 
ing radiation. Fluxes for the left and right beams are 

/out.L = Iin + QiuCOs{A{9k+K{9k))) + UinSm{‘i{6k + K{6k))) 
+ Vin5{9k) sm{2{9k + K{9k))), (A6) 

/out,R = lin-Qin COs{4:{9k +K{9k))) -Uin Sm{4:{9k + K{9k))) 
- VinS{9k) sm{2{9k + K{9k))). (A7) 


APPENDIX B: POLAROASTROMETRIC 
SIGNAL TRANSFORMATION TO ROTATED 
REFERENCE SYSTEM 

For practical purposes we frequently need to transform the 
polaroastrometric signal to a new reference system rotated 
relative to the initial system by some angle '0, measured 


counterclockwise. Taking into account the transformation 
of vectors and the Stokes parameters in the rotated refer¬ 
ence system, we derived these transformations for the com¬ 
ponents of the polaroastrometric signal as follows: 

Sq' ={Sq COS 24) + s* sin 24)) cos 4) 

-b {tq cos 24 > + tn sin 24 )) sin 4), 
tf = — (sj cos 24 > + Su sin 24 )) sin 4) 

+ {t* cos 24) + tn sin 24)) cos 4), 
s*f =(—Sg sin 24) + s* cos 24)) cos 4) 

+ {—tq sm 24 ) + t*n cos 24 )) sin'i/), 
tf = — (—Sq sin 24 ) + Sn cos 24 )) sin 4) 

+ {—t*q sin 24) -b tu cos 24)) cos 4), 

where values in new reference system are noted with primes. 
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For polaricentres we get 


where 


{aSq + bsu) cos Ip + {atq + bt^) sin-)/), 

— (asq + bsu) sin')/’ + {atq + bt„) cos tp, 
(csq + dsu) costp + (ctq + dt„) s'mip, 

— {csq + dSu) sintp + {ctq + dtu) costp, 


q cos 2i/) 

q cos 2-i/) + M sin 2tp ’ 
u sin 2%p 

q cos 2ip -\-u sin 2'tp ’ 
—q sin 2'tp 

—q cos 2'tp + u sin 2'tp ’ 
—wcos 2'tp 

—q cos 2'tp -\- u sin 2'tp 


(B2) 


(B3) 


As one can see, sometimes it is impossible to transform 
the polaricentres because the denominators in the last four 
expressions may be close to zero. In such cases one should 
operate with the polaroastrometric signal itself. 

It is worth noting that when the polaricentres q and 
u coincide, transformation of their coordinates is a simple 
rotation: 


s' = s cos 'tp + t sin 'tp, 
t' = —ssin'tp + t cos 'tp. 


(B4) 



